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Abstract 

Recently, the capture and storage of CO2 have attracted research interest as a strategy to reduce the global emissions of greenhouse 
gases. It is crucial to find suitable materials to achieve an efficient CO2 capture. Here we report our study of CO2 adsorption on 
boron-doped Cgo fullerene in the neutral state and in the le~-charged state. We use first principle density functional calculations to 
simulate the CO2 adsorption. The results show that CO2 can form weak interactions with the BC59 cage in its neutral state and the 
interactions can be enhanced significantly by introducing an extra electron to the system. 



Introduction 

The continuous dependence on fossil fuel combustion for the 
generation of energy has dramatically increased the atmos- 
pheric CO2 concentrations over the last century. Despite 
concerns for global climatic changes and many attempts to 
sustainably generate energy, fossil fuel combustion continues to 
be the main source of electricity while releasing 13 Gt of CO2 
[1] to the atmosphere each year. Therefore CO2 capture and 
storage (CCS) technology is a promising solution to reduce 
atmospheric CO2 emissions [2]. Solvent absorption that is based 
on amines is the most common technology for the capture of 
CO2. However this method is criticized for its very high energy 
consumption and operational limitations such as corrosion, slow 



uptake rates, foaming and large equipment. Hence there is a 
huge interest in solid adsorbent materials for CCS [3-6]. In past 
few years metal organic frameworks (MOFs) have emerged as 
solid CO2 adsorbent materials due to their tuneable chemical 
and physical properties. 

Particularly, there is growing interest for metal free carbon- 
based nanomaterials for gas adsorption. Carbon-based nanoma- 
terials such as fullerene, carbon nanotubes and graphene offer 
excellent thermal and chemical stability as CO2 adsorbents 
[7,8]. Heterofullerenes are fullerene structures in which one or 
more cage carbon atoms are substituted by heteroatoms [9]. In 
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addition to the properties mentioned above, which are inherent 
to carbon-based nanomaterials, heterofuUerenes also offer 
excellent tuneable chemical and physical properties [10]. Gas 
adsorption on heterofuUerenes is an appealing subject. B. Gao 
et al. [11] studied CO2 adsorption on calcium decorated Cgo 
fuUerene and F. Gao et al. [12] studied O2 adsorption on 
nitrogen-doped fullerene. 

Boron-doped C^o fullerenes are one of the most structurally 
stable heterofuUerenes [9]. Guo et al. synthesized B-doped C^q 
fullerenes for the first time, in microscopic amounts by laser 
vaporisation [13]. Zou et al. [14] demonstrated the synthesis of 
B-doped Cgo fullerene by using radio frequency plasma-assisted 
vapour deposition. Recently Dunk et al. [15] introduced a 
method to produce BC59 directly from exposing Cgo fullerene 
to boron vapour. Wang et al. [16] stated that substituting a 
single C atom of the C^q fullerene with a B atom does not cause 
a significant distortion in the cage structure. The net change in 
the dihedral angle due to the doping is only 1.6% and Kurita et 
al. [17] predicted that due to the similarity between the C-B 
bond and the C-C bond, the changes in the bond lengths are 
less than 5%. Therefore the BC59 fullerene has a similar struc- 
tural and thermal stability as C^o fullerene. Despite the 
numerous study results, which confirm the structural stability of 
B-doped €50 fullerene, very little studies have been done on 
applications of B-doped fullerene. Here, for the first time we 
report a study about the CO2 adsorption on B-doped Cgo 
fullerene, in which a single C atom is replaced with a B atom. 

Sun et al. [8] predicted an enhanced CO2 adsorption on \e~- 
and 2e~-charged boron nitride sheets and nanotubes, which 
show very little chemical affinity towards CO2 in their neutral 
state. Also Sun et al. [18] showed that chemical interactions 
between boron-carbon nanotubes (B2CNT) and CO2 can be 
enhanced by introducing extra electrons to the system. The 
enhanced interaction of CO2 with adsorbent materials by elec- 
tron injection has been further proved by Jiao et al. [19]. There- 
fore, we will investigate the CO2 adsorption on BC59 fullerene 
in both the neutral and the le~-charged states. 

Computational Details 

First-principles density functional theory (DFT) calculations 
were carried out to study CO2 adsorption on the BC59 cage. The 
BC59 structure was fully optimized in the given symmetry. The 
calculations were carried out at B3LYP [20-22] level of theory 
while using the split valance polarized basis set 6-31G(d). 
B97d [23,24] with the same basis set was used for calculations 
when non-covalent interactions are predominant. The CO2 
adsorption on BC59 was studied in the neutral state and in the 
le~-charged state. The electron distribution and transfer were 
analysed with MuUiken population analysis method [25]. 



The adsorption energies were calculated using the following 
equation. 

^ads = ^C02/BC59 " (^BC59 + ^€©2 ) ' (1) 

where E^^^ is the adsorption energy, Eqq^^q^^ is the total 
energy of the BC59 cage with a CO2 molecule adsorbed and 
^BC59 ^C02 energies of the isolated BC59 cage and 

CO2 molecule, respectively. For a favourable adsorption the 
calculated adsorption energy should have a negative value. To 
provide more accurate results for the chemi sorption energy the 
counterpoise corrected energy [26,27] was also calculated. 

The transition state was located by using the synchronous 
transit-guided quasi-Newton (STQN) method [28,29], which 
was then fully optimized by using the Berny algorithm at the 
B3LYP/6-31G(d) level. The optimized transition structure was 
used for IRC calculations at the same level of theory [30,31]. 
All calculations were carried out by using the Gaussian 09 
package [32]. The GaussView 5 package [33] was used to visu- 
alize the optimized molecular structures, molecular orbitals and 
charge distributions. 

Results and Discussion 

The substitution of a C atom in the C50 fullerene by a B atom 
causes a charge transfer between C and B atoms, which results 
in an unbalanced charge distribution in the fullerene cage. The 
unbalanced charge distribution forms B-C complex sites for the 
adsorption of CO2 (Figure 1). Here we considered two possible 
sites for the CO2 adsorption: the B-C atomic site between two 
hexagonal rings (HH B-C site) and two identical B-C sites 
between a hexagonal ring and pentagonal ring (HP B-C site). 

Adsorption of CO2 on uncharged BC59 
fullerenes 

According to our simulation results, the CO2 molecules can 
only form weak interactions with BC59 cage in its neutral 
state. The physisorption energy is a weak -2.04 kcal/mol 
(-4.1 kcal/mol for B97D/6-31G(d) calculations) and the weak 
interactions are mainly van der Waals interactions between the 
CO2 molecule and the adsorbent. The CO2 physisorbed con- 
figuration is shown in Figure 2. The CO2 molecule sits parallel 
to the boron-carbon plane of the BC59 fullerene cage. The B - Q 
and C - Q bond distances are 3.25 A and 3.71 A, respectively. 
The CO2 molecule undergoes very slight structural changes 
upon physisorption on the uncharged BC59 fullerene cage. The 
O-C-0 angle is slightly bent to 179.7° and the changes to the 
C=0 bond lengths are negligibly small. The doped fullerene 
cage hardly undergoes any structural change. The charge 
transfer between CO2 and BC59 is only 0.008e. 
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Figure 1 : Sites for CO2 adsorption on BC59. The B and C atoms of HH 
B-C and HP B-C sites are represented as 'ball and bond'-type and the 
rest of the atoms are represented as 'wireframes'. Atom colour code: 
grey, carbon; pink, boron. 




Figure 2: Configuration of physisorbed CO2 on neutral BC59. Atom 
colour code: grey, carbon; pink, boron; red, oxygen. 



Effects of charges on the structure 

Kim et al. [34] predicted that €596" should be a stable entity 
because of the isoelectronic configuration with C^q. This claim 
is further validated by experimental observations by Dunk et al. 
[15]. The Mulliken charge analysis and the electron density 
distributions of the lowest unoccupied molecular orbitals 



(LUMO) are adopted to assess the influence of changing the 
charge state of BC59. Figure 3 shows that the LUMO of the 
neutral BC59 is noticeably concentrated on the B atom and the 
neighbouring C atoms. Furthermore experimental results of Guo 
et al. [13] showed that boron doping creates an electron deffi- 
cient site at the B atom. This suggests that an additional elec- 
tron added to the system will be accepted by the B atom. This 
hypothesis is consistent with theoretical predictions of Kurita et 
al. [17] and Xie et al. [35], who stated that the doped B atom in 
Cgo fullerene acts as an electron acceptor. The comparison of 
the Mulliken population analysis of the neutral and the le~-state 
of BC59 proves that the negative charge introduced to the 
system is essentially accepted by the B atom. The Mulliken 
atomic charge of the B atom in the BC59 structure in the neutral 
state has changed from 0.138 to 0.012 upon the introduction of 
the negative charge, while as shown in Figure 4 the charges on 
the C atoms are not changed significantly. 




Figure 3: LUMO of neutral BC59. The orbitals are drawn at an isosur- 
face value of 0.02. The colours of the orbitals: red, positive wave func- 
tion; green, negative wave function. Atom colour code: pink, boron; 
grey, carbon. 



CO2 adsorption on BC59 fullerene in the 1e - 
state 

Next we studied the CO2 adsorption on a le~-charged BC59 
cage. The results confirm that the negatively charged BC59 
fullerene exhibits a stronger interaction with CO2. Unlike the 
neutral BC59, for which the interaction with CO2 molecule was 
only physical, here the charged BC59 forms a substantial chem- 
ical interaction with CO2 causing the molecule to undergo 
significant structural deformations. A stable CO2 adsorption is 
observed at the HH B-C site. The chemisorption energy of 
-15.41 kcal/mol (-64.48 kJ/mol) (-13.48 kcal/mol with BSSE 
correction) agrees well with the ideal range of chemisorption 
energy (40-80 kJ/mol) for a good CO2 adsorbent [36]. 
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Figure 4: Mulliken charge distribution of (a) neutral BC59 and (b) 1e -BC59. Tlie atoms are sliaded based on tine cliarge distribution on eacli atom. 
Tine comparison suggests tliat the most notable charge transfer is on the B atom. 



J 



The CO2 molecule undergoes considerable distortion upon 
chemically adsorbing on the le~-charged BC59 fullerene. A 
C=0 bond of the CO2 molecule is broken when one oxygen 
atom forms a bond with the boron atom (which will be referred 
as Oa in the following discussion and the other oxygen atom as 
Ob) and the C atom of the CO2 molecule forms a bond with the 
C atom on the HH B-C site of the cage structure. The linear 
0-C-O bond of CO2 is bent to 128.0° in the adsorbed form. 
The C=Ob bond which is originally 1.169 A (experimentally 
1.162 A [37]) is elongated to 1.208 A, while the length of the 
C-Oa bond is expanded to 1.336 A.The adsorption site of the 
BC59 fullerene also undergoes considerable stretching. The HH 
B-C site is protruded outwards by about 0.05 A. The B-C bond 
of the HH B-C site has stretched from 1.496 to 1.672 A. The 
Mulliken population analysis shows that a charge transfer of 
0.42 has occurred from the BC59 fullerene to the CO2 molecule. 



Comparison of the charge distribution on BC59~ before 
(Figure 4b) and after (Figure 5 c) CO2 adsorption, confirms that 
the injected electron is occupied by the CO2 molecule. 

The higher adsorption energy and the significant distortions in 
the structure confirm a stronger interaction between CO2 mole- 
cule and negatively charged BC59 than its neutral state. These 
interactions can be explained due to the Lewis acidity of CO2, 
which prefers to accept electrons [18]. On the other hand the B 
atom of the BC59 becomes less positively charged upon the 
addition of an extra electron. Therefore it becomes more likely 
to donate electrons to the CO2 molecule leading to stronger 
interactions between the two molecules. 

Figure 6 shows the minimum energy pathway for the adsorp- 
tion from the physisorbed state to the chemisorbed configur- 




Figure 5: (a) CO2 chemisorption and (b) transition structure for CO2 chemisorption on 1e -charged BC59. Atom colour code: grey, carbon; pink, 
boron; red, oxygen, (c) Charge distribution after CO2 chemisorption. 
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ation. We performed frequency calculations on the optimized 
transition structure, which confirms that it is a first order saddle 
point and hence an actual transition structure. From this figure, 
the activation barrier for the chemisorption is estimated to be 
13.25 kcal/mol (55.43 kJ/mol). The low barrier of the reaction 
indicates that the reaction is energetically favourable. 
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Figure 6: Intrinsic reaction pathway for CO2 cliemisorption on 
1e~-charged BC59 from tine pliysisorbed configuration. Tine total 
energy = 0 point corresponds to the total energy of Eco2 + Ebc59~^- 



For the desorption step, the removal of the added charge will 
decrease the stability of the bond between CO2 and the doped 
fullerene. The thermodynamic analysis of the reaction shows 
that the CO2 chemisorption is spontaneous only for tempera- 
tures less than 350 K. Therefore we suggest a method of manip- 
ulating the charge state and the temperature of the system for 
adsorbent recycling. Charging the system can be achieved by 
electrochemical methods, electrospray, and electron beam or 
gate voltage control methods [8]. 

Conclusion 

By using DFT calculations we have studied the adsorption 
mechanisms of CO2 on a C50 fullerene cage, in which a single 
C atom is substituted by a B atom. Our calculation results show 
that the BC59 cage, in its neutral state, shows a low chemical 
interaction with CO2 molecule, which only physisorbs with 
^ads ^ ~2.04 kcal/mol. However CO2 adsorption on the BC59 
can be significantly enhanced by injecting negative charges into 
the structure. The CO2 molecule chemisorbs on the le~-charged 
BC59 with ^ads ^ -15.41 kcal/mol. This study suggests that we 
can conclude le~-charged BC59 cage structure is a promising 
CO2 adsorbent. 
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